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Abstract 

We derive the gradient expansion for the diagonal elements of the Bloch propagator up to order 
f\^. We consider hamiltonians with position dependent effective mass and slowly varying local 
potential in spatial dimensions d = 1, 2, 3, 4. In our approach, we use the semiclassical h expansion 
of the density matrix in [Phys. Rev. A72, 022508 2005] with the help of Laplace transform to 
derive such gradient expansion in the one-body potential and also in the effective mass distribution. 
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I. INTRODUCTION 



Quantum mechanical system with a spatially varying effective mass have attracted a 
lot of attention and inspired intense research activities during recent years. They are for 
instance very useful in the study of the physical properties of various microstructures and 
semiconductor interfaces in condensed matter. Special applications of these models are 
carried out in the study of electronic properties of semiconductors [it, quantum wells and 
quantum dots. 2| Jl9| 3He clusters, jsl quantum liquids, 0] graded alloys and semiconductor 
heterostructures, [Sj . . . etc. These applications stimulated a lot of work in the literature on 
the development of methods and techniques for studying systems with mass that depends 
on position . On the other hand, the Bloch propagator or its Fourier transform namely 
the Green function are of prime significance since they contain all quantum mechanical 
informations on the system. To our knowledge the existing list devoted to the study on the 
Bloch propagator of quantum systems involvi ng p osition dependent effective mass is very 
short and is limited to the one dimensional case 10|, jsl. This has motivated the present work 
to obtain approximate analytical expression for the propagator or at least for its diagonal 
parts in spatial dimension d. The resulting expression may constitute a good test for the 
exact results. To this end, we make use of a semiclassical approximation, to derive simple 
expressions for the diagonal elements of the Bloch propagator when the effective mass is 
allowed to be position dependent. 

In a short, there exists a variety of procedures to generate the semiclassical h expan- 
sion with a final objective being the h expansion of the density matrix. We mention, for 
instance, the partition function method of Wigner-Kirkwood and further developped by 
Bhaduri and collaborators (see Ref. ll| and references cited therein)., the Kirzhnits expan- 
sion using commutator formalism [see Ref. and references cited therein] .i .H and 
the purely algebraic method introduced by Baraff and Borowich 14(] and further developed 
by Grammaticos and Voros 15|], based on the Wigner transform of operators. The latter 
method is particulary suitable for position dependent mass hamiltonians. These authors 
derived the semiclassical h expansions for the 3-dimensional one particle density and also 
for other densities of physical interest when kinetic energy operator of one-particle hamilto- 
nian comprises a spatially dependent effective mass. Later on, we have generalized [gI such 
expansion up to order /i^ for systems with effective mass distribution and reduced dimen- 
sionality, i.e, = 1,2 dimensions. It should be noted that corrections of order fi^ generate 
not only second order gradient corrections in the one-body potential and also second order 
gradient corrections in the effective mass distribution. Here, we are interested in obtaining 
the gradient expansion of the diagonal elements for the Bloch propagator m. d = 1,2 and 
4 spatial dimensions for hamiltonians with position dependent mass (the result in (i = 3 is 
already known ). 

The rest of the paper is organized as follows. In Section 2 we briefly recall some basic 
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definitions concerning the use of tlie Blocli propagator and its relationsliip to tlie density 
matrix. Starting from the semiclassical h expansion for the particle density, we derive in 
Section 3 the corresponding expansion up to second order in h for the diagonal elements of 
the Bloch propagator in d = 1, 2 and 4 spatial dimensions for hamiltonians with position 
dependent mass. A general analytical expression is found for the latter quantity in terms 
of the dimension d of the considered space. Section 4 provides an illustrative numerical 
example. Finally, a conclusion is given in Section 5. 



II. BASIC CONCEPTS 

Consider a system of N noninteracting fermions with spatially varying effective mass 
m*(T^) moving in a smooth potential U{1^). Throughout the present study, we shall be 
working with the one-body Hamiltonian given by 



H = -^'^.f{l^)^ + U{1^) (1) 

where /(^) = mo/m*^!^) denotes the ratio of the free particle mass uiq to the position 
dependent effective mass m*(l'") and where we use, as done in the majority of work on the 
subject, the symmetric ordering form of mass and momentum in the kinetic energy term of 
H. 

Let (pn{~>^) be the eigenf unctions of H and €n the corresponding eigenvalues, i.e; 
Hipn{l^) = En^Pni/f^), at zcro temperature the single-particle density matrix of the system 
7^'), is given by 

^0 = E^n(^)^"(^')0(A - e^) (2) 

n 

where A is the Fermi energy and B(x) is the Heaviside step function which allows to restrict 
the sum over occupied states only. 

Given the above density matrix, the Bloch propagator C("r^, "P"'; defined as [see for 
instance 



iii] 



C{1^, /3) := (T^l exp(-/3if) = 5^<(7^)^„(T^0 exp(-/3£,) (3) 

n 

can be obtained through the Laplace transform result 

C{'^,1^'][5) = (5 dAe~^V(^,^') (4) 
Jo 

It should be noted that, in quantum statistics and thermodynamics, /3 is an inverse of 
temperature: /3 = l/ksT with ks is the Boltzmann constant but if we now replace (3 in 
equation (3) by /3 — > it/h, the resulting propagator K^l^ ,l^']t) describes the propagation 
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of the single particle from 7^' — )■ ^ in time t. However in the subsequent analyis (3 is to be 
viewed as a complex parameter. The interest in the Bloch propagator is that it contains all 



quantum mechanical informations 21|, 11|, 22|, from which the density matrix p("^, in 
Eq. (2) may be obtained by suitable inverse Laplace transform, that is 

P(^, = ^ / e^'^^^^^f^d(3, c > (5) 

c—ioo 

where, to carry out the complex integration in equation (5), the parameter (3, as stated 
before, is considered as a mathematical variable which is taken to be complex. 
Putting "P"' = ^ in equation (4), we evidently get 

POO 

/3) = /3 / dXe-^^p{l^) (6) 







where henceforth C{1^;(3) denotes the diagonal elements of the Bloch propagator, called 
also the Slater sum, and p{l^) being the local particle density. 

In the above spirit, all the results are formally exact, we shall in the next section use 
them within the framework of gradient expansion. Notice that, for the case of constant 
effective mass hamiltonians the local version of equation (5) has been used to calculate 



the density from the Bloch propagator j20|. In the present context dealing with position 
dependent hamiltonians, we invert the procedure and we shall use equation (6) to get the 
Bloch propagator since, as stated before, the gradient expansions of the local density are 
known. 



III. GRADIENT EXPANSION FOR THE DIAGONAL ELEMENTS OF THE 
BLOCH PROPAGATOR FOR SPATIALLY VARYING EFFECTIVE MASS HAMIL- 
TONIANS IN DIMENSIONS D=l,2,3,4. 

In this section explicit h expansions will be presented for the diagonal elements of the 
Bloch propagator through the use of equation (6). For that, we directly use the h expansions 
of the particle density p("r^) derived inj^. In one spatial dimension it is given up to order 
h"^ by [see equation (A5) of reference j^] 
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Pd=i {x) 



2mn 



1 1 



3271 P \dx 



(A - V)-'^' + 



487r \ dx^ fdxdx ' 



1 



dV 



327r V c/x 
1 



{\-V) 



-5/2 



5(A-\/) + 



1 



247r \ (ix y 



dX 



(7) 



Here the potential V{^) is related to the one-body potential U{1^) in Eq. (1) by V"(7") = 
U{1^) + o^^^/("r^) and 5 is the Dirac distribution. Notice that, here we do not include 



the spin degeneracy (factor of two for spin half particles) in the expression of the particle 
density in (7) as was done inQ. Next, It is easy to verify that, the equation may be simplified 



to 



Pd=i 



2mn 



(A - l/)+i/2 ^ (A _ + 



1 1 



2mo I 167r ^dx 



d^V IdfdV 
2^^ + 



d {X - vy^^^ e {X - V) 

dX 

^d{X-Vy^^^9{X-V) 
dX 



247r \ dx"^ f dx dx 



+ 



1 rdv 



24:71 \ dx 



d^x-vy^^^e{x-v) 



(8) 



which can be alternatively rewritten as 

Pd=i {x) 



2mo 
vr V W 



{X-Vf'^9{X-V) + 



f9f 



1 1 



dl 

2mo I 167r \dx 



2^^ + ■' 



'd{X-Vf^^e{X-V) 



dX 



1271 \ dx"^ f dx dx 



' d^X - Vf'^^ e {X - V) 



+ 



1 fdV 



127r V dx 



'd\x-vY^^^e{x-v) 



dX^ 



(9) 
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The reason that we have displayed the result in the form (9) is that, for the case of a constant 
effective mass our expression of the density reduces exactly to the one given by equation (12) 
of jol in one spatial dimension, where the latter density was obtained throught a different 
semiclassical method namely the Kirzhnits expansion. 

To obtain the local Bloch density, we use for the density the form in (8) rather than (9) 
and substitute it into Eq. (6) and then we use the following property of Laplace transforms 



|l7| applied to each term of the expansion 



/ dXe^^^[{\-v{-^)Y e{x-v{-^))\ 

Jo 



(10) 



and using obvious notations, we then obtain 

1/2 



C^d=i(^;/3) 



1 df dV 
f dx dx 



1 + 



h^f \3 fidf 



24mo \ 4: \f dx 



(3+ 



- 2 



d^V 

dx"^ 



dy_ 

dx 



(11) 



Let us now come to the two dimensional case and we write down the corresponding h 
expansions of the local density given by Eq. (28) of ^ 



Pd=2 



:{\-v) e(A- V) + ^ 



2nh?f 



487r 



f 



f 



(12) 



Plugging Eq.(12) into (6) with the use of Laplace transform identity (10), we find for the 
local Bloch propagator 
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C^d=2(^;/3) 



27r fffPj 



24mo 



/ 



f 



For the d = 3 case, the expression of the local density up to order h'^ [see for instance [16| 
and references cited therein] 

3/2 



Pd=3 




{X - Vf^ Q{X - V)+ 



f 



{x-vf^e{x-v)+ 



(A - vy^^^ e(A - V) -i (^^vy (A - vy''^^ e(A - v) 

(14) 



6 



Substituting Eq. (14) into (6) followed by the use of (10), we get 



3/2 



\2TTh?f(3 

f 



1 + 



V 



24mo 
f3^ + (^Vyf3 



7 (ff^ 
4 / 



f 



/3+ 



(15) 



It is interesting to observe that the results given respectively in equations (11), (13) and 
(15) can be written down in terms of the dimensionality d of the space as follows 



d/2 

^0 \ ^-I3V{1^) ^ 



(d-2) 



1 + 



24mo 



[d-l] 



/3 + 



/ 



V 



(3^ + {^VYf3' 



(16) 



We have looked whether the above equation is valid for higher dimensions or at least for 



d = 4. For that, we have make use of the semi classical approach in Ref. 15|] to write down 
the local density up to order h'^ in d = 4 dimension. We have obtained 



=4(7*) 



f 



V7 

/ 



i\-v)ei\-v) + 



mo 



48712^2/ 



f 



(17) 



Here ^ stands for the gradient in four dimensions. Having the density, we follow the same 
derivations as done for d = 1,2,3 to obtain the corresponding local Bloch propagator, it 
yields 



X 



1 + < 3 

24mn 



/ 



/ 



/3 + 



Vf.VV 

f 



V 



{3' + 
(18) 



It is easy to check that, setting c? = 4 in Eq.(16), one recovers the result in (18). Hence our 
expression given in (16) holds true for d = 1,2, 3, 4. Notice that, in equation (16) unlike the 
position dependent mass terms, the remaining terms involving gradients of the potential do 
not depend on the dimension d of the considered space. 
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The interested reader may have noticed that the local Bloch density in Eq. (16) is 
expressed in terms of V^(l^) = U(l^) + -^^"^ fi^) and for practical use it is interesting to 
re-express it in terms of the original one body potential U{1^). Upon substitution, equation 
(16) becomes then 



d/2 

"mo \ ^-i3U(i^) 



2TTh^fl3 



id-2) 



V 



d/2 

mo \ ^-i3U(-^) 



2Trh?f/3 



id-2) 



f 



1 + 



1 + 



24mo 



(3 + 



V 



+ i^V)'f3' 



_/3 '^-^2/(1^) 

8mn •' ^ ' 



24mo 



(d-l)^ + 3 (Yf 
f 



[l-d) 



f 



/3 + 



V 



+ i^V)'f3' 



(19) 



where in getting the second form, use has been made of the Taylor expansion up to order /i^ 
of exp(— /3g^ V 2/) ~ 1— /Sg^^^/- Since, the terms involving and ^'^V are of order 
h"^, we need only replace them in Eq. (19) by their leading terms, i.e 
and U + Oih?) . This leads finally, after simple rearrangements, up to order h'^ 



Cdi^;P) 




id + 2) 



f 



/3 + 



(20) 



where we recall that /("P") = mo/m*(7'"). As can be seen the above local Bloch density 
receives explicit contributions from the spatially varying effective mass ■m*i'r^) not only 
at zero order, through the term proportional to /""'^^ and also from terms of order 
proportional to and ^2 

/. The above equation is the main result of the present study. 
For d = 3, our expression reduces to the one obtained long time ago in the context of nuclear 
physics [23 1 . 



IV. NUMERICAL EXAMPLE 

In this section, we want to test numerically the importance of the position depen- 
dent effective mass terms in the resultant propagator (20). For that, we need as an 
input a given effective mass m*and potential U. As an illustration we focus on one- 
dimensional systems with mass distribution m*(x) and we choose f/(x) so that equation 



(1) posseses exact analytical solution. In Ref. |2J], Alhaidari solved exactly equation 
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(1), for a large class of potentials, by means of an elegant method called point canonical 

transformations (PCT). Let us briefly sketch this technique. Under the following (PCT), 

_ i_ _ 1 

y = J {f{x)) 2 dx and (pn{x) = (/(x)) * ipniv), the Schrodinger equation (1) with spa- 
tially mass distribution m*{x) = mo//(x) and potential U{x), is mapped to a Schrodinger 
equation with a constant mass mo, so that. 



+ u{y) My) 



U{x) = U{y) + 



8m(x) 



1 d?m(x) 



7 f dm{x) 
im?{x) \ dx 



and Er, 



En^niy), with 
Taking for the con- 



m{x) dx^ 

stant mass problem the harmonic oscillator potential U{y) = mou'^y'^ /2, Alhaidari obtained 
the solutions for the Oscillator class for a given m*(x). Let us now take the specific mass 
distribution used in \24\. 



m{x) = niQ 



7 + 



x^ 



m(±oo) = mo 



(27) 



from which we get f{x) = ((1 + a;^) + with 7 being a real constant parameter, 
one then obtains ?/ = x + (7 — l)arctan{x). Note that when 7 = 1.0 equation (27) gives a 
constant effective mass m(a;) = mo- In terms of f{x) and its derivatives, the above potential 
U (x) reads then 



Uix) 



rriQijj 



[x + (7 — l)arctan(x)]' 



d'fjx) 
dx"^ 



+ 



/(^) 



df{x) 
dx 



(28) 



In Figure 1., we plot the ratio /(x) = mQ/m{x) from Eq.(27) as a function of 
the spatial coordinate x for values of 7 = 0.6, 0.8, 1.0. Setting the parameters cu 
and /3 to 1, we plot in Figure 2. the Bloch density Cd=i{x; (3) of Eq. (20) as a 
function of x for values of 7 = 0.6, 0.8, 1.0. Here we consider particle of unit free 
mass mo = 1, and the Planck constant is set to 1. As is shown in this figure the 
presence of spatially varying effective mass may lead to important deviations locally 
with respect to the constant mass case. We also display in Figure 3. the second 
order h"^ term of the Bloch density, which using Eq. (20) is given by 6Cd=i{x; (3) = 



1/2 



-I3U 



4 \f dx I f dx^ 



24mo \2nh^fl3 J 

Deviations from the constant effective mass case are clearly exhibited at this order. 



/3 + 



1£_ dU 
f dx' dx 



9 I d?U_ 
^ ' dx^ 



V. CONCLUSION 

In the present study we have derived the gradient expansion up to order W'^ for the 
diagonal elements of the Bloch propagator for hamiltonians with position dependent mass. 
Our main result is summarized by equation (20) valid for spatial dimensions d = 1,2,3,4. 
This result is valid for arbitrary one-body potential U{1^) and spatially varying effective 
mass m*{^). 
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FIG. 1: A plot of position dependent effective mass ratio /(x) = mo/m*{x) = ((l + x^) /(7 + a^^))^ 
.The solid curve is at value 7 = 1, the dot curve is at value 7 = 0.8 and the dash curve is at value 
7 = 0.6 
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FIG. 2: A Plot of The local Bloch density Cd=i{x; (5) with a; = 1, /3 = 1. Solid line correspond to 
7 = 1, doted line 7 = 0.8 and dashed line 7 = 0.6 . 
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FIG. 3: The second order contribution to the Bloch density in figure 2 6Cd=i{x; (3) for values of 
7 = 0.6,0.8, 1.0 corresponding respectively to dashed, doted and solid lines. 



We have explicitly shown that, for d = 1 and constant effective mass the expressions of 
the semiclassical densities derived respectively within the algebraic method due to Baraff 
and Borowich ij] and the Kirzhnits expansion are identical up to order h?. This would be 



certainly true for higher dimensions. 

The results we have obtained would constitute useful approximation to the exact cal- 
culation of propagator or Green functions for hamiltonians with spatially varying 
mass. 

An interesting extension of the present study will be the derivation of the semiclassical 
expansion for the non-diagonal elements C(7", 7"'; /3) of the Bloch propagator. 
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